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Abstract 

We develop a class of new kinetic data structures for colli- 
sion detection between moving convex polytopes; the per- 
formance of these structures is sensitive to the separation 
of the polytopes during their motion. For two convex poly- 
gons in the plane, let D be the maximum diameter of the 
polygons, and let s be the minimum distance between them 
during their motion. Our separation certificate changes 
0(log(D/s)) times when the relative motion of the two poly- 
gons is a translation along a straight line or convex curve, 
0{\/D/s) for translation along an algebraic trajectory, and 
0(D/s) for algebraic rigid motion (translation and rotation). 
Each certificate update is performed in 0(log(D/s)) time. 
Variants of these data structures are also shown that exhibit 
hysteresis — after a separation certificate fails, the new cer- 
tificate cannot fail again until the objects have moved by 
some constant fraction of their current separation. We can 
then bound the number of events by the combinatorial size 
of a certain cover of the motion path by balls. 

1 Introduction 

Collision detection is an algorithmic problem arising in 
all areas of computer science dealing with the simula- 
tion of physical objects in motion. Examples include 
motion planning in robotics, virtual reality animations, 
computer-aided design and manufacturing, and com- 
puter games. Often the problem is broken up into two 
parts, the so-called broad phase, in which we identify the 
pairs of objects we need to consider for possible collision, 
and the narrow phase in which we track the occurrence 
of collisions between a specific pair of objects [14]. (In 
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the spatial database literature, these are also called the 
filtering and refinement phases, respectively [22].) For 
the broad phase, almost all authors use some kind of 
simple bounding volumes for the objects themselves, or 
for portions of their trajectories in space or space-time, 
so as to quickly eliminate from consideration pairs of 
objects that cannot possibly collide. The narrow phase 
is more specialized, according to the types of objects 
being considered. 

The simplest objects to consider are convex polytopes 
(polygons in the plane, or polyhedra in 3-space), and 
this case has been extensively considered in the litera- 
ture [17, 18, 20, 10, 5]. More complex objects are then 
broken up into convex pieces, which are tested pairwise. 
Algorithmically, the convex polytope intersection prob- 
lem is a special case of linear programming; in two and 
three dimensions even more efficient techniques have 
been developed in computational geometry, that can be 
applied after a suitable preprocessing of the two poly- 
topes [6, 7]. The methods, however, that have proven 
to work best in practice exploit the temporal coherence 
of the motion to avoid doing an ab initio intersection 
test at each time step. Not surprisingly, the collision 
detection problem is closely related to the distance com- 
putation problem. Since the distance between two con- 
tinuously moving polytopes also changes continuously, 
many well-known collision detection algorithms, such as 
those of Lin and Canny [16, 17], Mirtich [18, 19, 20], and 
Gilbert et al. [10] (see also [4]), are based upon tracking 
the closest pair of features of the polytopes during their 
motion (which, of course, implies knowledge of the dis- 
tance between the polytopes). The efficiency of these 
algorithms is based on the fact that, in a small time 
step, the closest pair of features will not change, or will 
change to some nearby features on the polytopes. 

Though it is hard to imagine how one can do bet- 
ter than tracking the closest pair of feature when the 
polytopes are in close proximity, such tracking seems to 
be unnecessarily complicated when the polytopes start 
moving further from each other. Indeed most of the 
above authors suggest performing first a simple bound- 
ing volume (box or sphere) test on the two polytopes, 
and only if that fails entering the closest feature pair 
tracking mode. In this paper we consider a number 
of general techniques that allow us to perform colli- 
sion detection between two moving convex polytopes 
in a way that is sensitive to the separation between the 
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polytopcs. In order to properly quantify the separation- 
sensitivity of our methods, we view the collision detec- 
tion problem in the context of kinetic data structures 
(or KDSs for short), introduced in [2, 12]. 

In the kinetic setting we assume that the instanta- 
neous motion laws for our polytopes are known, though 
they can be changed at will by appropriately notify- 
ing the KDS. Our sampling of time is not fixed, but is 
determined by the failure of certain conditions, called 
certificates. In our case these are separation certificates, 
which prove that the two polytopes do not intersect. 
The failure of a separation certificate need not mean 
that a collision has occurred; it can simply mean that 
that certificate has to be replaced by one or more oth- 
ers, still proving the non-intersection of the polytopes. 
A good KDS is compact if it requires little space, re- 
sponsive if it can be updated quickly after a certificate 
failure, local if it adjusts easily to changes in the motion 
plans of the objects, and efficient if the total number 
of events is small. Our kinetic collision-detection data 
structures have all these properties; they maintain only 
a small constant number of certificates, and the cost for 
processing a certificate failure or a motion plan update 
is at most polylogarithmic (in the combinatorial size of 
the polytopes). 

A number of the papers referenced above make the 
claim that their algorithms are efficient because "if the 
sampling interval is small enough, then the cost of up- 
dating the closest pair of features is 0(1)." This is a 
difficult statement to attach a fully rigorous meaning 
to — exactly how small the time step must be to guar- 
antee this condition depends a great deal on both the 
polytopes and the speed and complexity of their mo- 
tion. In the kinetic setting we can give the notion of 
efficiency a more satisfactory theoretical definition, by 
focusing on the maximum number of certificate failures 
we may have to process for polytopes and motions of a 
certain complexity, rather than on the adequacy of any 
absolute unit of time. 

The key contribution of this paper is to develop a 
class of new kinetic data structures for collision detec- 
tion between convex polytopes, where the efficiency of 
the structure can be analyzed in terms of natural at- 
tributes of the motion. Given two moving convex poly- 
gons in the plane, let /i = min{n, y/D/a}, where n is 
the combinatorial complexity of the polygons, D is their 
maximum diameter, and a is their minimum separation 
during the entire motion. In Section 5 we develop a 
KDS where the number of events (certificate failures) is 
O (log fi) when the relative motion of the two polygons 
is translation along a convex trajectory (for example, 
a straight line), O(n) for translation along a algebraic 
trajectory, and 0(/j, 2 ) for algebraic rigid motion (trans- 
lation and rotation) . Thus we see how the nature of the 
motion, as well as the proximity of the two polygons, 



affect the complexity of the collision detection problem. 

In contrast to this, the closest pair of features of two 
polygons can change fi(n 2 ) times under an algebraic 
rigid motion, no matter what their separation is. For 
'intermediate separation' situations, when the bound- 
ing boxes of two objects intersect, but their distance is 
still Q(D), our methods will perform much better than 
other extant methods for collision detection. The per- 
formance of our methods interpolates smoothly those of 
the bounding box and closest pair of features techniques 
mentioned above, as the separation varies. In this in- 
termediate distance range our methods are also directly 
useful for non-convex objects, as such objects can be 
bounded by their convex hulls. 

We attain these distance-sensitive bounds by con- 
structing certain novel outer approximating hierarchies 
for our polytopes, whose structure is of independent in- 
terest. These hierarchies provide a series of combina- 
torially simpler and simpler shells, as we move away 
from the polytope. For two polytopes in proximity the 
hierarchies locally refine so as to provide a separation 
certificate. 

Hidden in the above 'O' bounds are factors depend- 
ing on the algebraic degree of the motion. Again, when 
the polytopes are in close proximity, it is clear why a 
'wiggly' motion should be more costly than a smooth 
one. But why should it be so when they are further 
away? This has motivated us to develop structures that 
exhibit hysteresis — where, after a certificate failure has 
occurred, no other certificate failure can happen until 
the objects have moved by some constant fraction of 
their current separation. Using these structures, we are 
able to bound the number of events by the combinato- 
rial size of a certain cover of the motion path by balls 
(Section 6), somewhat reminiscent of [21]. 

We believe that the KDSs shown here are of theo- 
retical and practical interest. A basic tool for all our 
structures are certain outer approximation hierarchies 
for convex polytopes and their Minkowski sums — a topic 
which we believe to be of independent interest (Sec- 
tion 3). The distance- and motion-sensitive bounds we 
give are novel and, to our knowledge, the first such to 
be ever presented. 

It was surprising to us that even for the simple set- 
ting of two moving convex polygons, there is much that 
is novel and interesting to say; in fact, many challeng- 
ing open questions remain. Though our exposition is fo- 
cussed on the two-dimensional case, we do not anticipate 
significant obstacles in extending our results to three di- 
mensions; we briefly describe some preliminary results 
in Section 7. We expect that our kinetic structures will 
lead to improved practical algorithms for convex shapes, 
and we have already started an implementation of our 
algorithms. In a companion paper [3], we discuss a dif- 
ferent set of kinetic collision techniques applicable to 
non-convex shapes. 
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2 Models of motion 

In our model, each object is a closed rigid convex poly- 
gon, whose motion is described by a moving orthogo- 
nal reference frame: a point o(t) and two orthogonal 
unit vectors x(t),y(t), whose coordinates are continu- 
ous algebraic functions of t. (Such moving frames do 
exist, and are flexible enough to approximate any mo- 
tion to any order and accuracy, for a limited time. Note 
however that an algebraic rotation is necessarily of non- 
uniform angular velocity, and can cover only 0{d) full 
turns, where d is the degree of the entries.) Each vertex 
of the object is assumed to have constant algebraic co- 
ordinates (X, Y) relative to the frame (o, x, y); so that 
its position at time t is o(t) + Xx(t) + Yy(t), also an 
algebraic function of t. 

As we shall see, the certificates used in each of our 
KDSs have the general form F(t) = G(pi(t), . . . ,p n (t)). 
Here pi(t), . . . ,p m (t) are the positions of m vertices (m 
will be a small constant), possibly on different objects; 
and G is some algebraic function. Then F itself is an 
algebraic function, which means we can compute the 
next time t when F(t) — 0, exactly, and compare any 
two times, at finite cost, within an appropriate arith- 
metic model. Moreover, the number of zeros of F in 
any finite interval is bounded by its algebraic degree, 
which is the degree of G times the maximum algebraic 
degree d among the motion coordinates. So, for, ex- 
ample, the same vertex triplet cannot become collinear 
more than 2d times. 

In conclusion, if the complexity (algebraic degree and 
coefficient size) of the motions is bounded, each certifi- 
cate can fail at most 0(1) times during a single motion, 
and the cost of computing and comparing the failure 
times is also 0(1). We will use these last two facts, 
and only these two facts, throughout our analyses. Our 
results do not require any other combinatorial prop- 
erties of algebraic motion — for example, that an alge- 
braic path can be decomposed into a constant number 
of convex sub-paths, as in the companion paper [3]— 
and therefore apply to a wider range of pseudo-algebraic 
motions [2, 12]. 

3 Polygon Approximation Hierarchies 

Our collision detection algorithms are based on outer 
approximation hierarchies for the convex polygons in- 
volved, or for their Minkowski sum. All these hierarchies 
tile the space exterior to the polygon so as to simplify 
the combinatorial structure of the polygon as one goes 
away from the polygon. A well-known example of such 
a hierarchy in computational geometry is the Dobkin- 
Kirkpatrick hierarchy [7] . This hierarchy is not by itself 
adequate for our purposes, however, as it is not sensi- 
tive to the distance away from the polygon; any approx- 
imation in the hierarchy can have vertices arbitrarily far 




Figure 1. A boomerang hierarchy. 

from the original polygon. We will retain a key property 
of this hierarchy, namely the fact that that each 'coars- 
ening step' is performed by removing an existing edge of 
the current approximation and extending outwards its 
two neighboring edges till they meet. However, in order 
to guarantee that we get the distance-sensitivity that we 
require, we will enrich the set of 'available edges' by in- 
troducing an additional set of degenerate edges around 
the polygon, initially all of zero length. In general, the 
number of these degenerate edges will be proportional 
to the number of original edges in the polygon. 

Though one normally thinks of these approximating 
hierarchies as being constructed from the boundary of 
the polygon towards the outside, it is actually advan- 
tageous to visualize this process in reverse. Let P be 
our convex polygon of n edges. We start from the out- 
ermost approximation Pq, which for our purposes is al- 
ways a (not necessarily axis-aligned) bounding rectan- 
gle of P. The space between P and Pq comprises up to 
four non-convex polygons, each consisting of a concave 
chain and two additional edges. Following Hershberger 
and Suri [13], we call these polygons boomerangs. The 
last coarsening step, if undone, corresponds to 'cutting 
off a corner' of Pq through the reintroduction of an- 
other (possibly degenerate) edge of P. This step splits 
one of the top-level boomerangs into a triangle and two 
smaller boomerangs. This process of cutting off corners 
is continued recursively, until all edges of P have been 
reintroduced. See Figure 1. Structurally, the hierar- 
chy consists of four binary trees, where each internal 
node corresponds to a triangle and each subtree to a 
boomerang. In our hierarchies, each of these trees will 
have depth O(logn). A variety of outer approximations 
for P can be defined by removing a subtree of triangles 
rooted at each of the top-level boomerangs. 

Before we discuss various choices of degenerate edges, 
let us develop some notation and terminology. The 
apex of a boomerang is the unique vertex not on the 
boomerang's concave chain. The height of a boomerang 
is the distance from its apex to its concave chain; this is 
the maximum distance from any point in the boomerang 
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to the chain. The level of either a boomerang or a trian- 
gle is its depth from the root in the appropriate binary 
tree; there are at most 4 • 2 l boomerangs at level i. Fi- 
nally, let Pi denote the i-th envelope of P, defined as 
the union of P and all level-i boomerangs; Pj is itself a 
convex polygon surrounding P. 

Let D denote the diameter of P. The following lemma 
implies that the envelopes in any boomerang hierarchy 
of P are reasonably close to P. We omit the easy proof 
from this abstract. 

Lemma 3.1. In any single level of any boomerang hi- 
erarchy there are 0(^J D/ s) boomerangs of height at 
least s. 

A line £ that does not intersect P can intersect at most 
one of the four top-level boomerangs in any boomerang 
hierarchy of P. Moreover, if £ intersects a boomerang, 
then it intersects at most one of its two sub-boomerangs. 
In fact, these two observations hold for any convex curve 
(bending away from P). These simple observations es- 
tablish the following useful lemma. 

Lemma 3.2. Any convex curve disjoint from P inter- 
sects at most one triangle in each level of any boomerang 
hierarchy of P. 

We observe that the triangles in a boomerang hierar- 
chy for P tile the space between Pq and P. (It is not 
hard to extend this tiling to be the complement of P 
in the plane by allowing a few infinite triangles. This 
is almost identical to the construction of a binary space 
partition tree [9].) For a point x outside P, the trian- 
gle in this tiling that contains x provides us with useful 
information about the position of x with respect to P: 
the base of the triangle is a side of P separating x from 
P, while the height of the associated boomerang is an 
upper bound on the distance from x to P. 

3.1 The Compass Hierarchy 

For the compass hierarchy we introduce 0(n) zero- 
length edges into P, whose outer normals form a regular 
recursive lattice on the unit circle, in the standard com- 
pass directions (E, N, W, S, NE, NW, SW, SE, etc.). 
In the top |~log 2 ri] levels of the compass hierarchy, each 
boomerang is subdivided into two smaller boomerangs 
and an isosceles triangle. In the remaining levels, if 
any, each boomerang is subdivided by a line through its 
median edge, as in a standard Dobkin-Kirkpatrick hi- 
erarchy. The resulting hierarchy has at most 2 |~log 2 ra] 
levels. 

Lemma 3.3. (a) A level-i boomerang in the compass 
hierarchy of P has height 0(D/2 i ). 

(b) The compass hierarchy of P contains 
Oy^/D/s log(£)/s)J boomerangs with height at 
least s. 




Figure 2. The Dudley construction. Solid vectors introduce 
degenerate edges; dashed vectors introduce degenerate vertices. 

(c) Any convex curve at distance s from P intersects 
0(log(D / s)) triangles in the compass hierarchy of P. 

(d) If two polygons P and Q are distance s apart, 
their compass hierarchies contain disjoint approxima- 
tions P and Q, each with 0(log(D / s)) edges. 

Klosowski et al. [15] define the "fc-DOP" or discrete 
orientation polytope of an object to be the bounding 
polytope whose facets are normal to a fixed set of k 
'compass' directions. Klosowski et al. construct a hier- 
archy of bounding volumes for any object by computing 
the object's fc-DOP, decomposing the object into a con- 
stant number of pieces, and recursively constructing a 
hierarchy for each piece, using the same value of k at all 
levels. (See [1, 11, 24] for similar bounding volume hi- 
erarchies.) In contrast, the compass hierarchy consists 
of a nested sequence of fc-DOPs with k = 4, 8, 16, ... , 
all for the same object. 

3.2 The Dudley Hierarchy 

Our second boomerang hierarchy is based on a result 
of Dudley [8] on approximating convex bodies in arbi- 
trary dimensions by polytopes with few facets; hence, 
we call it the Dudley hierarchy. Let S be a set of n reg- 
ularly spaced points on a circle of radius 2D, centered 
inside P. For each point x £ S, let n(x) be its nearest 
neighbor on P. If n(x) is a vertex of P, we introduce 
a zero-length edge at n{x) whose outer normal vector 
is n(x) — x. Otherwise, we introduce a new degenerate 
vertex at n(x) whose external angle is zero. See Fig- 
ure 2. We then create a boomerang hierarchy, starting 
with the bounding box Pq, by recursively subdividing 
each boomerang by a line through its median (possibly 
degenerate) edge (and possibly through other collinear 
edges). The resulting Dudley hierarchy has depth at 
most |~log 2 (2n)] . 

Lemma 3.4. (a) A level-i boomerang in the Dudley 
hierarchy of P has height 0{DI^) [8]. 
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(b) The Dudley hierarchy of P contains 0(y/D/s) 
boomerangs with height at least s. 

(c) Any convex curve at distance s from P intersects 
0(\og(D/ s)) triangles in the Dudley hierarchy of P. 

(d) If two polygons P and Q are distance s apart, 
their Dudley hierarchies contain disjoint approximations 
P and Q, each with 0(\og(D/s)) edges. 

4 Mixed Hierarchies 

The hierarchies introduced so far are tilings of the free 
space around a single convex polygon. Since we are in- 
terested in systems with two (or more) moving convex 
polygons, we would like to extend some of these notions 
to hierarchical tilings of the free part of the configu- 
ration space generated by the motion of two polygons. 
Our general plan is to associate (explicitly or implic- 
itly) a particular separation proof with each tile of such 
a tiling. As the two polygons move and their configu- 
ration crosses out of its current tile, some certificates 
will fail and a new separation certificate will have to be 
generated. 

Let P and Q be our two moving polygons. If P and 
Q are only translating with respect to each other, then 
the configuration space remains two-dimensional. It is 
well-known that in this case the free space is the com- 
plement of the convex polygon Pffi (— Q), the Minkowski 
sum of P with the negative of Q. Since again our two- 
dimensional configuration space is the exterior of a con- 
vex polygon, then all the hierarchies presented above 
can be directly used, if we arc willing to construct this 
polygon. Such a direct approach, however, is less at- 
tractive when rotation is allowed, as then the configura- 
tion space can have quadratic combinatorial complexity. 
Also, for applications where we have multiple moving 
polygons (though we do not discuss techniques for this 
many-body problem in this paper), we want to 'recy- 
cle' as much as possible pieces of the hierarchies built 
for each of the polygons, rather than having to build a 
separate hierarchy for each pair. 

We describe below one such hierarchy for tiling the 
complement of the Minkowski sum of two convex poly- 
gons into triangles and parallelograms, which we call 
the mixed hierarchy. (Without loss of generality, we 
focus on describing this hierarchy for P © Q, as op- 
posed to P © (-Q).) It is motivated by the theory of 
mixed volumes of convex bodies [23] , and it has the ad- 
vantage that it changes in a very regular way as the 
polygons rotate. This makes it possible to maintain a 
non-intersection certificate for two moving convex poly- 
gons by separating the motion into a translational part, 
corresponding to a point moving around in the plane, 
and a rotational part, corresponding to a deformation 
of the triangle or parallelogram containing the point. 
At certain discrete events this tile and a neighboring 
tile are deleted and replaced by two other tiles cover- 




Figure 3. Inserting an edge into the mixed hierarchy. 

ing the same area, much like a Delaunay flip. Thus, 
we can maintain a separation proof by maintaining one 
tile, or a small working set of tiles, around the current 
configuration point. 

Given two boomerang hierarchies for P and Q, the 
mixed hierarchy is defined by starting from Po © Qo 
and then interleaving the 'corner cutting' operations 
leading to P and Q. Different mixed hierarchies may 
be obtained for different interleavings of these corner 
cutting operations. The main difference from a stan- 
dard boomerang hierarchy arises because the corner of 
P or Q being cut next may no longer be a corner at 
all in the Minkowski sum hierarchy: between two con- 
secutive sides of P we can have many sides of Q, and 
vice versa. To be concrete, say the next cut is to add 
a side e of P, but its neighboring sides in the P hier- 
archy are now separated by a chain of Q edges in the 
mixed hierarchy built so far. Note that these neighbor- 
ing sides must already be present in the mixed hierarchy, 
as the joint corner cutting sequence is consistent with 
that for P. To insert e into the mixed hierarchy, we 
fist find the place where e fits, according to slope, in 
the chain of Q-edges that replaced its corner. We parti- 
tion that Q-chain at that point, and then translate the 
two subchains inwards, as illustrated in Figure 3. The 
two chains come to rest when they encounter the points 
where e meets its two neighboring edges in P. Thus this 
process adds to the mixed hierarchy exactly the same 
triangle that it added in the P hierarchy, as well as sev- 
eral parallelograms based on the edges of the Q chain. 
After all the corners have been cut, the space between 
Po © Qo an d P © Q is tiled with one copy of each of 
the original triangles used in the P and Q hierarchies, 
as well as several 'mixed' parallelograms. 

Figure 4 shows a full mixed hierarchy for two small 
convex polygons. The corner cutting order is indicated 
by the numbers next to the edges of the polygons. 

As we mentioned earlier, a nice aspect of the mixed 
hierarchy is its behavior under polygon rotation. The 
rotation changes the interleaving of the edges of P and 
Q around their Minkowski sum. For example, an edge 
of P can pass in slope ordering one of the edges of the 
Q-chain that it splits. The primary effect of this change 
on the mixed hierarchy is a simple change, which we 
call a a page turn, akin to a Delaunay flip: a triangle 
and a parallelogram exchange positions. Additionally, 
at some finer level of the hierarchy, one parallelogram 
may either collapse to a line segment and be destroyed, 
or a new parallelogram may be created and expand from 
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Figure 4. A mixed hierarchy for two convex polygons. 

a line segment. See Figure 5. 

Lemma 4.1. If the corner cutting operations for poly- 
gons P and Q of size n and m, where to < n, are inter- 
leaved according the level of the boomerang destroyed 
at each step (in the P or Q hierarchy respectively) , then 
the size of the mixed hierarchy is 0((m + n) log to). 

Lemma 4.2. If P is stationary and Q makes a single 
full rotation, then the number of page turns that will 
happen is 0(mnlogm). 

Lemma 4.3. Given a cell c of the mixed hierarchy, a 
point p £ dc, and some auxiliary structures of linear 
size for P and Q, the neighboring cell of c containing 
p can be computed in O(logn) time. The same applies 
for the cells covering c, when cell c is destroyed during 
a page turn event. 

If the boomerang hierarchies of P and Q satisfy dis- 
tance properties such as those in Section 3, then similar 
results will hold for the tiles of the mixed hierarchy. 
For example, if we mix two compass or Dudley hierar- 
chies, we can reduce the time bound in Lemma 4.3 from 
O(logn) to 0(\og(D/s)), where s is the distance fromp 
to Pffi Q. 

In general, however, the mixed hierarchy does not sat- 
isfy the equivalent of Lemma 3.2. Although any line 
disjoint from P © Q hits only one triangle per level, it 
may hit a large number of that triangle's neighboring 
parallelograms. Since the compass hierarchy uses the 
same evenly-spaced cutting directions for every poly- 
gon, at most one parallelogram appears next to any tri- 
angle in the mixed compass hierarchy, so this problem 
is avoided. Unfortunately, we have no similar guarantee 




Figure 5. Events in the mixed hierarchy under rotation. 

for the Dudley hierarchy, so our bound in that case is 
much weaker. 

Lemma 4.4. (a) The mixed compass hierarchy 
and the mixed Dudley hierarchy of P © Q have 
0(^D/s log(D/s)) and 0(y/D/s) vertices with 
height at least s, respectively. 

(b) Any convex curve at distance s from P®Q inter- 
sects 0(\og(D/s)) and 0(yjD/s) triangles in the mixed 
compass hierarchy and mixed Dudley hierarchy of P®Q, 
respectively. 

5 Separation-Sensitive Data Structures 

Several authors have proposed algorithms to maintain 
the closest pair of features between two convex polygons 
[10, 17, 18]; these algorithms can easily be transformed 
into a kinetic data structure with constant update time. 
There are at least two alternative approaches that lead 
to the same performance. We could maintain an in- 
ner common tangent between P and Q, i.e., a line that 
touches the boundaries of both objects, but separates 
their interiors. Alternately, we could maintain a sepa- 
rating edge e of one polygon, along with the vertex v of 
the other polygon closest to the line through e. Unfortu- 
nately, for all three approaches, the worst-case number 
of kinetic events is quite high — 0(n) if the polygons are 
only translating, and 0(n 2 ) if they are also allowed to 
rotate. Moreover, these lower bounds can be achieved 
while the polygons are arbitrarily far apart. 

In this section, we describe several new kinetic data 
structures that maintain a separation certificate be- 
tween two moving convex polygons, where the cost and 
number of events depends on the distance between the 
polygons. Our data structures are loosely based on the 
algorithm of Dobkin et al. [6] for detecting intersections 
between preprocessed convex polygons or polyhedra. 

Let us first establish some notation. P and Q are 
convex n-gons with diameter at most D. We let s denote 
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the current separation (geometric distance) between P 
and Q at a given time, and a the minimum of s over 
the entire history of the motion. Finally, we let fi = 
min{n, y^D/a}. 

5.1 One Point, One Polygon 

We illustrate our basic approach by considering the 
special case where P consists of a single point p. 
We construct either a compass or Dudley hierarchy 
Qoi Qi, Q2, ■ ■ ■ for Q. Each triangle at level i in this 
hierarchy has an inner edge, which is an edge of Qi+i, 
and two outer edges, which are subsets of edges of Qi. 

At any moment, we maintain the active triangle A 
containing the point p. There are two types of certificate 
failures; see Figure 6(a). If p crosses one of the outer 
edges of A, we can identify the triangle it enters in 
constant time. On the other hand, if p crosses the inner 
edge of A, it could either collide with Q or pass into a 
triangle at some deeper level of the hierarchy. We can 
check for actual collision in O(l) time be seeing if p lies 
on the line segment eC\Q. Otherwise, we search one level 
at a time for the new active triangle A'. If A is at level 
i and A' is at level j, we find A' in 0(j —i) = O(logn) 
steps. Alternately, if we use binary search, we can find 
A' in time 0(log(logn — i)). 

If the point is moving along a convex path (curved 
away from Q), then by Lemma 3.2 at most one triangle 
in each level is ever active. For the Dudley and compass 
hierarchies, the level of the deepest active triangle, and 
thus the number of active triangles, is 0(\ogp). We 
easily observe that the total time spent updating the 
active triangle is also only 0(log/i). 

Now suppose the point is moving along some other 
algebraic path. For the compass hierarchy, the number 
of events is 0(p log p) , and for the Dudley hierarchy, the 
number of events is 0{p). Both upper bounds follow 
directly from the number of triangles that can contain 
a point at distance a or higher from Q (Lemmas 3.3(b) 
and 3.4(b)). As in the case of convex motion, these are 
also upper bounds on the total update time. 

Theorem 5.1. As p moves algebraically about Q, we 
can maintain a separation certificate maintained in 
0(min{loglogrt,log(Z)/s)}) time per event, using 0(n) 
space and preprocessing time. Both the number of 
events and the total update time is 0(log)i) if p moves 
along a convex curve and 0{n) otherwise. 

Despite its good performance, the KDS just described 
is somewhat wasteful. Only the inner edge of the active 
triangle separates the moving point from the polygon, 
so why worry about the other two edges? This obser- 
vation motivates the following lazy variant of our KDS, 
which has exactly the same performance bounds as The- 
orem 5.1 in the worst case, but is likely to be more effi- 
cient in practice. 




(a) (b) 

Figure 6. Events for a point moving around a polygon, (a) Main- 
taining the active triangle, (b) Maintaining a separating edge. 

Instead of a triangle, we maintain a single separating 
edge ei of some envelope Qi. We update the separating 
edge only when p passes through the line containing 
ej. Again, there are two types of events: if p hits the 
edge ei, we have a stab event, and otherwise, we have a 
push event. See Figure 6(b). Stab events are handled 
exactly as in the previous structure: first check for a 
real collision, and if no collision has occurred, find a 
new separating edge at a deeper level in the hierarchy. 
After a push event, one of the edges of Qi adjacent to ej, 
say e\, is now a separating edge. However, the structure 
of the hierarchy ensures that e\ is actually a subset of 
an edge ej of some coarser envelope Qf, we take e 3 to 
be the new separating edge. 

5.2 Two Translating Polygons 

Now consider the case of two convex polygons P and Q 
which are translating along algebraic paths. It suffices 
to consider the case where Q is fixed and only P moves. 
Detecting collisions between P and Q is equivalent to 
detecting collisions between a single moving point and 
the static Minkowski sum P (— Q), so the bounds 
in Theorem 5.1 immediately applies to the case of two 
translating polygons. 

This structure is unsatisfactory, however, since it re- 
quires us to construct a decomposition of the Minkowski 
sum P (B (— Q). If we want to collisions among several 
convex n-gons using this approach, we need 0(n) space 
for every active pair of polygons. In this section, we 
describe a modification of our previous data structures 
that use a separate hierarchy for each polygon, so that 
we only need 0(n) space for each polygon, plus constant 
space for every active pair. 

Our first approach is to maintain the active cell c of 
the mixed hierarchy of P (— Q) containing the cur- 
rent configuration point. Whenever the configuration 
point leaves c, we can compute the new active cell c' 
in 0(\ogp) time (Lemma 4.3). We emphasize that it 
is not necessary to construct the entire mixed hierar- 
chy explicitly, but only separate hierarchies for the two 
polygons. If we build the mixed hierarchy out of the 
Dudley hierarchies of P and Q, convex translation can 
now cause 0(/i) events. If we use compass hierarchies 
instead, convex translation causes only 0(log u) events, 
but the event bound for more general translations goes 
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(a) 




00 

Figure 7. Events for two polygons undergoing (a) translation 
and (b) rotation. 

up to 0(fj,logfi). As in the one-point case, the event 
bounds also bound the total update time. 

Theorem 5.2. As P and Q translate algebraically, 
we can maintain a separation certificate maintained in 
0(log(min{n, D/s}) time per event, using 0(n) space 
and preprocessing time per polygon plus 0(1) extra 
space for the pair. Both the number of events and the 
total update time is 0(\og /i) if P moves along a convex 
curve relative to Q and O(n) otherwise. 

We can also achieve these bounds with a suit- 
able modification of our lazy KDS. Let Po,P±,... 
be a boomerang hierarchy for P, and Qo,Qi,--- a 
boomerang hierarchy for Q. As the polygons move, 
we maintain a separation certificate (ej,Uj), where 
is an edge of Qi and Vi is the closest vertex of Pi to 
&i (the line through a), or vice versa. We always use 
features of the same level in both hierarchies. Given a 
valid separation certificate (ej,Uj), we can compute ei- 
ther a finer certificate (ej+i, t>j+i) or a coarser certificate 
(ei_i,«i_i) in constant time, if one exists, by checking 
local neighborhoods of ej and Uj. 

The separation certificate expires when v\ crosses < e~t . 
As before, there are two types of events; see Figure 7(a). 
(In the following description, we will assume without 
loss of generality that m € Pi-) 

To handle a stab event, where m hits ej, we refine 
both envelopes one level at a time until they are disjoint. 
(Since a single refinement may introduce a zero-length 
edge at i>j, we may have to refine by more than one 
level.) If Pj and Qj are the coarsest disjoint envelopes, 
then either the edge of Qj containing ej or an edge of Pj 
containing Vi is a new separating edge. If «j is actually 
a vertex of P and it hits the edge of Q containing ej , we 
report a collision. 

After a push event, where Vi does not hit the edge 
ej, either an edge of Qi adjacent to e, or an edge of Pi 
adjacent to Vi is a separating edge between Pj and Qi 
(or possibly both), depending on which has the higher 
slope. After updating the separation certificate, we 
coarsen the envelopes one level at a time, computing 
a new separation certificate at each new level, until the 
next coarser envelopes intersect. 



Motion 


Compass 


Dudley 


One point, one polygon 


convex 
general 


OQogft) 
0(/ilog /z) 


0(\ogn) 
0(jm) 


Two polygons 


convex translation 
general translation 
rigid motion 


OQogn) 
0(/ilog n) 
0( M 2 logV) 


0{n) 
O(p) 



Table 1. The number of events and total update times of our 
KDSs, for different types of objects, motions, and boomerang 
hierarchies. Here, fj, = min{n, -J D/cr}, where D is the objects' 
maximum diameter and a is their minimum separation. 

Since we can refine or coarsen by one level in constant 
time, the total cost of either event is 0(log/i). The 
number of events is the same as for the mixed-hierarchy 
structure. 

5.3 Rigid Motion 

Now suppose P and Q are also allowed to rotate. As we 
mentioned earlier, the mixed hierarchy changes as the 
polygons rotate. If the active cell in the mixed hierarchy 
disappears due to a page turn, we can construct the new 
active cell in 0(log/i) time. Since we only maintain the 
active cell, page turns elsewhere in the mixed hierarchy 
cost us nothing. We do not have to worry about collapse 
events, since the configuration point will be 'squeezed 
out' of the cell before it finishes collapsing. 

Our lazy two-hierarchy data structure can now en- 
counter an additional event, called a roll, when one of 
the edges adjacent to becomes parallel to e,*; see Fig- 
ure 7(b). To handle a roll event, we keep 6j as the sep- 
arating edge, but now its nearest neighbor in the other 
envelope is a vertex adjacent to m. After we update the 
separation certificate, we then coarsen the envelopes as 
much as possible in 0(log/i) time, just as for a push 
event. 

For both approaches, we have the following bounds, 
if we use the Dudley hierarchy; the event bound for the 
compass hierarchy is slightly weaker. 

Theorem 5.3. As P and Q undergo algebraic rigid 
motion, we can maintain a separation certificate in 
0(log(min{n, D/s}) time per event, using 0(n) space 
and preprocessing time per polygon, plus O(l) space 
for the pair. Both the number of events and the total 
update time are 0(/u 2 ). 

5.4 Summary 

Table 1 summarizes the event bounds for our kinetic 
data structures for various types of objects, classes of 
motion, and boomerang hierarchies. 
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Figure 8. An inflated hierarchy. Only the outermost level of 
(dashed) connecting edges is shown. (Compare with Figure 1.) 

6 Inflation, Hysteresis, and Path- 
Sensitivity 

By further modifying the lazy variants of the kinetic 
data structures described in the previous section, we 
can obtain data structures that exhibit hysteresis: after 
any event, the configuration must change by a certain 
amount before the next event occurs. Hysteresis allows 
us to derive upper bounds on the number of events based 
on geometric properties of the path that, unlike our ear- 
lier results, do not depend on smoothness or algebraic- 
ity. 

For any convex polygon P and real number e > 0, the 
(outer) offset polygon P[e] is a convex polygon obtained 
by moving each edge of P outwards by a distance of 
e and moving the vertices along their angle bisectors. 
That is, each edge of P[s] is parallel to an edge of P 
and vice versa, and the distance between their two lines 
is e. For any point p on P[e], we have e < d(p, P) < 
ej sin(f?/2), where 6 is the minimum internal angle at 
any vertex of P. Zero-length edges in P induce positive- 
length edges in P[e]. 

As in the previous section, we first consider the spe- 
cial case of a single point p moving around a polygon Q. 
Let Qo, Qi, ... be the Dudley hierarchy of Q, and let 
£j = aD/2\ the approximation error of Qi given by 
Lemma 3.3(a). (We can obtain similar results using a 
Dudley hierarchy instead.) We define an inflated hier- 
archies Q' ,Q[, ■ ■ ■ , where Q[ = Qi[ei]. See Figure 8 
(but ignore the dashed edges for now). Since Qo is a 
rectangle, no envelope Qi has an acute vertex angle, so 
d{Q'i,Q) < (1 + V2)ei. 

As the point moves, we maintain a separating edge e 
of some uninflated envelope Qi. The separation certifi- 
cate expires when p crosses the line 4 e~ > . To compute the 
new separation certificate, we find the index j such that 
p lies outside Q'j but inside Q'j-i- For the new sepa- 
rating edge, we take the edge of Qj parallel to the edge 
of Q'j that intersects e . To keep the required space 
down to 0(n), we cannot precompute the breakpoints; 
instead, we compute each breakpoint e n Q' on the 



fly in 0(j) time when we need it. The total cost of an 
event is 0(log fi), and the total number of events is the 
same as for the uninflated compass hierarchy. 

Alternately, we can connect successive levels in the 
inflated hierarchy, as shown by the dashed edges in Fig- 
ure 8, decomposing Q' Q \ Q into a complex of triangles 
and trapezoids. Each cell in this complex has one inner 
edge and two or three outer edges. After any event, we 
locate the cell c in the inflated complex containing p; if 
the inner edge of c is an edge of Q'j, we take the cor- 
responding edge of Qj to be the new separating edge. 
The update time and number of events is the same as 
above. 

Lemma 6.1. After any event, if d{P,Q) — il(D/n), 
then P must move at least d(P, Q)/ (3 before the next 
event, where (3 = 2(1 + \/2) « 4.8284. 

The only time we do not obtain hysteresis is when 
the new separating edge is an edge of the actual poly- 
gon Q, or equivalently, when the point p lies inside the 
polygonal annulus Q[e]\Q for some e = O(Dfn). 

For any real k > 1, say that a circle C is n-clear if 
its radius is at most 1/k times the distance from its 
center to the polygon Q. A n-clear decomposition of a 
path 7r is a decomposition of tt into contiguous segments 
7i"i , 7T2, . . . , 7Tfc, such that each segment 7Tj is contained in 
a K-clear disk. The size of such a decomposition is the 
number of segments. 

Theorem 6.2. If p moves along a path tt whose min- 
imum distance to Q is £l(D/n), then the number of 
events is at most the size of the smallest n-clear decom- 
position of it, where n = 2/3 + 1 = 5 + 4^2 w 10.6569. 

The constants (3 and n are function of the minimum 
external angle of any envelope and the ratio £j/£,_i. By 
using more complicated bounding polygons as the out- 
ermost level of the hierarchy, and by letting the inflation 
offset grow more slowly, we can decrease (3 arbitrarily 
close to 1 and k arbitrarily close to 3. Somewhat para- 
doxically, however, these modifications increase both 
the update time per event and the worst-case number 
of events. In particular, using an inflated Dudley hier- 
archy instead of an inflated compass hierarchy doubles 
the value of (3, even though it leads to asymptotically 
fewer events in the worst case. 

For the case of two translating polygons P and Q, 
we construct separate inflated hierarchies for both P 
and Q and use a technique similar to the lazy structure 
in Section 5.2. The separation certificate consists of an 
vertex v of Pi and an edge e of Qi, or vice versa, for 
some i. When v hits e, we compute a new separation 
certificate based on which inflated envelope Q'j contains 
the vertex v. The resulting structure exhibits hysteresis 
and path-sensitivity (with different constants (3 and k) 
and still satisfies Theorem 5.2. 
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Finally, a further modification gives us rotational hys- 
teresis as well. That is, after any event, P must ei- 
ther move a distance of fl(d(P, Q)) or rotate by an an- 
gle of Cl(d(P,Q)) relative to Q before the next event 
occurs. This modification requires a lower bound on 
the external angles of any envelope, so we can use only 
the compass hierarchy, not the Dudley hierarchy. Ro- 
tational hysteresis implies that the number of events 
is bounded by the size of a K-clear decomposition of 
the path that the polygons traverse through the three- 
dimensional configuration space, for some constant k. 

We omit further details from this version of the paper. 

7 Conclusions and Open Problems 

As we mentioned in the introduction, we do no foresee 
any major obstacles to generalizing the two-dimensional 
results in this paper to three dimensions. We already 
have a few preliminary results, which we plan to develop 
further in the full paper. Generalizations of both the 
compass hierarchy and the Dudley hierarchy are easy to 
define, and at least the Dudley hierarchy provides good 
approximation bounds. After constructing the Dudley 
hierarchy of two polyhedra, we can maintain a separa- 
tion certificate in 0(polylog(D/s)) time per event, using 
an approach similar to the lazy structure in Section 5. 
The number of events is 0(/j, 2 ) for algebraic translation, 
and 0(/i 4 ) for algebraic rigid motion. 

Returning to our two-dimensional results, for vari- 
ous measures of kinetic efficiency, we have KDSs with 
good performance under that measure, but only at the 
expense of some other kinetic quality. It would be desir- 
able to find a single KDS that uses only linear space per 
polygon, has good event bounds under all three classes 
of motion (convex translation, algebraic translation, and 
algebraic rigid motion), and exhibits both translational 
and rotational hysteresis. 

Finally, we intend to work on extending the two body 
methods (narrow phase) of this paper to multiple mov- 
ing convex polytopes. This will require a kinetic struc- 
ture that implements the broad phase of collision detec- 
tion and determines which pairs of objects need to be 
passed on to the narrow phase. 

Acknowledgments: We wish to thank Julien Basch 
for fruitful discussions, and Pankaj Agarwal for making 
us aware of Dudley's result [8] . 
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